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Abstract. A brief summary of strangeness mile stones is followed by a chemical 
non-equilibrium statistical hadronization analysis of strangeness results at SPS 
and RHIC. Strange particle production in AA interactions at ^sjvjv > 8.6 GeV 
can be understood consistently as originating from the deconfined quark-gluon 
plasma in a sudden hadronization process. Onset of QGP formation as function 
of energy is placed in the beam energy interval 10-30A GeV/c. Strangeness 
anomalies at LHC are described. 
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1. Introduction 

We consider, in a chemical analysis of yields, the strange hadron production in A-A 
collisions at the SPS and RHIC. We find that these results can be well understood and 
consistently described within the Fermi type, chemical near-equilibrium, statistical 
hadronization model G] . We find that the deviation from equilibrium conditions 
are most pronounced at RHIC (^snn = 130 GeV, RHIC130) HHH|, with the conditions 
of hadronization strongly deviating from chemical equilibrium. The present analysis 
and reanalysis of the results obtained in the energy range 19.4 < ^snn < 8.75 GeV by 
NA49 OHIBlElUniHIlHl and WA97/NA57 Q31HI1E3E0 experiments is confirming 
these features characteristic of the formation of a new state of matter at SPS. 

The strongest difference between SPS and RHIC collision results is the order of 
magnitude enhancement in strange particle yield per participant, this result is not 
new but has found so far little attention 0] E] . The strangeness yield rise increases 
compared to SPS [5| Ej in a manner which is more spectacular than the increase 
in the total hadron multiplicity. This is an expression of the increase of the excess 
of strangeness (greater chemical nonequilibrium) at RHIC as compared to SPS. We 
establish this increase in strangeness yield with precision, which requires use of an 
appropriate analysis tool of hadron multiplicities, which allows reliable and precise 
description of the unobserved particle phase space in domains not yet experimentally 
accessible. 

At RHIC 130, each baryon participant leads to formation of nearly 9 strange quark 
pairs in Au-Au reactions. We show in section|3that this high yield is reached gradually 
as function of the available fireball thermal energy. This is confirming that this high 
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yield of strangeness is result of the same physical process in all reaction energies 
considered here, beginning with ^/snn = 8.75 GeV (corresponding to a 40A GeV 
beam on laboratory stationary target). We also present in section how the diverse 
properties of the hadronic fireball change as function of energy. The introduction of 
chemical nonequilibrium leads to a rather low hadronization temperature which we 
interpret in sectional where we discuss in depth at which energy domain the onset of 
QGP may be occurring. 

Before we turn to this part of our report a brief introduction to 'strangeness' 
the important mile stones of 10 years of experimental research are presented in the 
following section [5] This is followed in section [3] by a discussion of the mechanism of 
QCD based production of strangeness in a quark-gluon plasma (QGP). In particular 
we show that the gluon fusion processes in QGP is capable to populate within the 
reaction time the strange quark phase space ^I^IEHGOI, U P to, and even above 
the chemical equilibrium yield in the QGP phase |2(J| . 

The other QCD based process, light quark-antiquark fusion into strange quark 
pair, has been earlier shown to be too slow for strangeness equilibration in the 
QGP Thus the discovery of the dominance of the gluon fusion process has 

been the key stepping stone in establishing strangeness as an observable of QGP. The 
computation presented in section is assuming early thermalization, and following 
rapid collective evolution. This scenario is today viewed as a very likely reaction 
sequence, seen the behavior of the V2 flow parameter, especially so that it applies also 
to strange hadrons [2"3| . 

A very important aspect of the study of strange particles produced by QGP is the 
final state hadron production process. We present the essential elements of statistical 
hadronization in section These principles are employed in section in analysis of 
the RHIC and SPS results, and as noted a synthesis of our analysis and their impact 
for the search of the onset of QGP formation is offered in section^ In the final section 
0we argue that there is a very interesting future for strangeness in the LHC energy 
domain. 

2. Strangeness — a popular QGP diagnostic tool 

There are several different strange particles in nature allowing us to study several 
complementary physics questions. Denoting here light flavors with q = u,d we consider 
the stable particles: 

<f>(ss), K(qs), K(qs), A(qqs), A(qqs), S(qss), S(ps), Sl(sss), fi(sss). 

Moreover, the study of strange hadron resonances such as K*, £*, A* is a 
very promising avenue of forthcoming research, leading to direct determination of 
hadronization conditions. 

Many strange hadrons are subject to a self analyzing decay within several 
centimeters from the point of production, as is illustrated in figure^for the cascading 
S~-decay. Thus one tracking device, typically today a TPC combined with silicon 
vertex detector allows the measurement of the many particle yields of interest. The 
accessibility of strange hadron distributions has lead, in recent years, to an explosion 
of both experimental and theoretical interest. Moreover, production rates of hadronic 
particles and hence statistical significance even for low intensity beam is usually high. 

Multistrange hadrons can be formed in hadronization of QGP by 'cross talk' 
between quarks made in disjoint microscopic reactions |24|. An illustration of this 
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Figure 1. Illustration of the S — -production by an incoming, e.g., SPS Pb-beam, 
on a heavy laboratory target, and its subsequent decay. Dashed line depicts the 
invisible neutral A, emerging from the decay kink E~ — > 7r~A, and ending in the 
decay 'V of A — * p + 7r~ . 




Figure 2. Illustration of the cross-talk two step mechanism of strange hadron 
formation from QGP: inserts show gluon fusion into strangeness, followed by QGP 
recombinant hadronization. 

process is presented in figure [21 the inserts show gluon fusion processes gg — > ss |17|. 
which establish ample supply of strange quark pairs in the early stages of QGP 
evolution. In the ensuing hadronization process, quark recombination leads to 
emergence of particles such as here shown R(sqq), £l(sss), which otherwise could 
only very rarely be produced, considering that several rarely occurring processes 
would have to coincide. Enhancement of such particles is for this reason indicative of 
deconfincmcnt |17M25j . This mechanism leads to the expectation that the enhancement 
of strange antibaryons is progressing with strangeness content, as was recognized and 
discussed qualitatively already 20 years ago, see e.g., section 4 in |26| . 

Since the production of strangeness occurs predominantly in thermal gluon fusion 
gg — > ss, see section [31 the overabundant presence of strangeness in high energy 
relativistic reactions can be on theoretical grounds linked to presence of thermal 
gluons, and thus to QGP. Beyond the dominant processes gg — ► ss there is also, 
at 10-15% of the total rate, a contribution from qq — ► ss. This rate alone would 
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however not suffice to equilibrate strangeness in rapidly disintegrating QGP phase 
formed in relativistic heavy ion collisions |21M22| . 

The remarkable coincidence of scales involving the strange quark mass m s , and 
the critical temperature of hadronic phase transition T c , 

m s ~ T c which implies r s ~ tqgpj (1) 

allows the gg — ► ss production process to act as a clock for the collision reaction. 
Another interesting feature is that, in a baryon rich environment, we also have 
s > q |24| . which applies when the quark chemical potential is greater than the 
mass of the strange quark. In this condition a strange antibaryon enhancement is 
expected, comparing to the nonstrange antibaryons, including the anomaly A/p > 1. 
This should be the case at the lower energy limit of SPS. 

In this line of thought, even though there is a near u,d,s flavor symmetry at 
RHIC, we also expect (anti)hyperon dominance of (anti)baryons in this condition. 
Here, the mechanism is slightly different |20| : each baryon or antibaryon contains 
three quarks or, respectively, antiquarks, and it is hard to find in the flavor symmetric 
pool of quarks (and antiquarks) three non-strange quarks (or antiquarks), in fact the 
statistical probability for assembly of three nonstrange quarks is 1/3. Thus 2/3 of all 
produced baryons, or respectively antibaryons, will in limit of flavor symmetry carry 
strangeness, and this feature of course remains true in the future LHC environment. 

Multistrange anti-hyperon abundance ratios were studied extensively and we show 
in figure the three ratios: 

ss d _ IF + 0.5 S* (1530) ssd E~ + 0.5 S*(1530) 

isd ~ A + SO + 0.92^(1385)' sud ~ A + E° + 0.92 S*(1385) ' 

ssd + ssd _ + 0.5 1 * (1530) + Br + 0.5 S*(1530) 

sud + sud ~ A + EO + 0.92 S*(1385) + A + S° + 0.92 S*(1385) ' 

Above, we have explicitly identified the main resonance contributions, accompanied 
by the appropriate branching ratio. Given their large spin-isospin quantum numbers, 
these resonances despite their greater mass are as important contributors to the final 
yields, as is the ground state. For the A, we also note the presence of inseparable 
E° — > A + 7 contribution. 

It is practically impossible to separate decay contributions which dilute the 
intrinsic ratios ssd/ sud, ssd/ sud down to about 0.23 for the case of RHIC |31j. 
However, we see a clear enhancement (factor 3) over the p—p, p—p background for all 
A-A interactions here considered. Among these, the RHIC result is of most profound 
consequence as we shall address below, as it implies large excess of strangeness yield 
(compare figure IT4"Jl . The enhancement of the antibaryon rations seen in the p-A 
ratio (middle field of figure |3J) is not seen for baryons (see right field of figure [3Jl . A 
simulated p-A enhancement can arise from the annihilation of antihyperons produced 
by the reacting matter in the in surrounding spectator nuclear matter. Namely, the 
annihilation cross section for singly strange hyperons is notably larger than for double 
strange hyperons. In consequence the A are more depleted than 3. 

To avoid the semblance of p-A enhancement process it is important to consider 
how the yield of individual hadrons behaves compared to expectations based on a 
cascade of N-N interactions. This view on the strange hadron enhancement has been 
presented by the CERN-SPS WA97 collaboration for 158A GeV Pb collisions with 
laboratory Pb target, and these results have been confirmed and extended in the more 
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Figure 3. Ratios of hyperon abundances produced in high energy interactions. 
Recent experimental results from E3EH1EH1E1E1EIE21E] ■ 



recent NA57 measurements. In figure^] we show the rise of the enhancement with the 
strangeness content, with p-Be reaction system used as reference. These results show 
yields per participant, where a participant is understood to be an inelastically reacting 
'wounded' nucleon. These results confirm the QGP prediction of enhancement growth 
with a) strangeness b) antiquark content. 

Before these studies of strangeness enhancement were carried out within the 
central rapidity region by NA57 and WA97 experiments, it has been shown that 
central rapidity is indeed the kinematic domain where the new physics is occurring. 
The NA35II experiment obtained A yield for 20(M GeV Sulfur projectile on laboratory 
Sulfur target As can be seen in figure [SJ the antibaryon yields arc localized 

in the central CM (center of momentum) rapidity region. The background p-p 
production result, shown by squares is obtained scaling up the yield by the number 
of nucleon participants. We see that the central rapidity is in fact the domain of the 
large enhancement as would be expected in a reaction picture in which a fireball of 
hot matter is formed. 

Strangeness yield as function of rapidity is obtained evaluating (s + s) — 
1.6 A + 1.6 A + 4Kg, and is seen in figure El Again, the measured result is reflected 
at CM rapidity and the squares provide a hadron multiplicity scaled p-p result as 
basis. Thus, the enhancement, seen in figure is expressing how much faster than 
nonstrange hadron strangeness increases in S-S interactions compared to p-p. 

One would expect naively, in such a comparison, a reduction rather than an 
enhancement, since in a cascade of conventional hadronic reactions, it would be easier 
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Figure 4. Abundance enhancement with respect to the yield in p— Be 
collisions, scaled with the number of 'wounded' nucleons, results of the WA97 
collaboration 1131 . 




y 

Figure 5. A, in NA35II 200/1 GeV S-S interactions 1321 . Background (squares) 
from multiplicity scaled N-N reactions forward (open circles) points: reflection at 
the CM rapidity. 



to produce less massive pions, than strangeness, remembering that cascading hadrons 
degrade in energy. There can be no doubt seeing this result that even in the relatively 
small S-S reaction system there is a new physical mechanism of strangeness production 
at work. In absolute number, the yield of strangeness, seen in figure^! is rather large, 
above 25 ss pairs. Our predictions that this could be happening should QGP phase 
be formed were classified 20 years ago at the LBL Heavy Ions Studies among 'exotica'. 
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Figure 6. CM-rapidity reflected 1.6A + 1.6A + 4K S 200A GeV S-S multiplicity 
scaled 200A GeV p-p results drawn from NA35II experiment 1331 . 
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Figure 7. Central rapidity spectra of single strange Kg, A and A reported by 
the experiment WA85 I34j. Antibaryon to baryon ratio _Ra = A/A is independent 
of mj_. 



Aside of particle yields, an important further key piece of physical evidence is 
the shape of transverse mass m± = J m 2 + p 2 ^ spectra. The experimental results 
cannot be understood within a conventional hadron cascade reaction picture. When 
the (central rapidity) m± spectral distributions are fitted to an exponential shape, one 
finds that the inverse m± slopes show a rather precise baryon-antibaryon universality, 
even though at CERN-SPS significant baryon number asymmetry prevails. This was 
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Figure 8. Central rapidity spectra of strange baryons A(sud), E.~(ssd), Q(sss) 
and their antibaryons A, S, Q reported by the experiment WA97 I35I . 

first discovered in S-induced collisions by WA85 experiment [32], see figure This is 
confirmed by the WA97 experiment in Pb-Pb Interactions |35]> see figure with the 
slopes given in tabled The inverse slopes of the baryon and antibaryon spectra are to 
a very great precision (1%) the same. Also, within the error, the slopes of A and 5 are 
the same. Results of the experiment NA57 indicate that the more precisely measured 
fl and spectra are also yielding the same T± . We understand the kaon-hyperon slope 
difference to be result of the collective explosive flow and the different range of 
m± considered evaluating the m± slopes. 

These mj_-slope results, found in different reaction systems, cannot be accidental, 
also since a different pattern of behavior is clearly seen in p-p and p-A reactions. The 
identity of baryon and antibaryon m± slopes observed in A-A reactions prove that 
strange baryons and antibaryons are produced by the same mechanism in nuclear 
collisions, and do not suffer rescattcring in a dense baryon-rich medium after their 
production. A systematic experimental [5J3 and theoretical study of the m± 
spectra for the central rapidity Pb-Pb reactions at projectile energy 158A GeV 
as function of reaction centrality confirms that Kaons and strange hyperons and 
antihyperons are born simultaneously and do not undergo rescattering. To the best 
of our understanding the only physical mechanism which can lead to this result is a 
QGP fireball breakup into free streaming hadrons, as would be expected in sudden 
hadronization of a QGP fireball [551 1591 13U1 HT] . 

In figure[51 we show a survey of the SPS-158A GeV transverse slopes. The inverse 
slope parameter T is related to the intrinsic thermal parameters T t [ in the source and 
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Table 1. The inverse slope Tj_ of mj_ spectra seen in figure^] Pb— Pb interactions 
at 158A GeV, WA97 experiment ESI- 
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Figure 9. Inverse slopes of particle spectra obtained for Pb-Pb interactions at 
projectile energy 158A GeV. 



local flow velocity v t {, for m,px S> T by the usual Doppler formula 



T^i±^T tf ^Ji±^T tf . (2) 

Vi - v t 2 f V l - «tf 

3. Kinetic description of strangeness production in QGP 

We view a QGP fireball as consisting of quarks and gluons constrained by the external 
'frozen color vacuum'. Within the deconfined domain these particles can move and 
undergo collisions, and reactions, within a volume much larger than the usual nucleon 
size. The reactions occur at energies characterized by thermal equilibrium and thus 
the kinetic process like this is often called 'thermal production' in order to distinguish 
these soft reactions from hard kinetic processes occurring prior to thermalization of 
the incoming particle momenta. These first collision processes play a decisive role in 
providing abundance of charm, and heavier flavors, in quark-gluon plasma. 

However, 'thermal production' can mean something entirely different in other 
related work. Namely, this term has also been used to describe the statistical yield 
of particles in chemical equilibrium. Both usages coincide only when the dynamics of 
the reaction process allows for the kinetic processes to establish chemical equilibrium. 
This, in general, will not be the case in heavy ion collisions at relativistic energies. 
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The generic angle averaged cross sections for (heavy) flavor s, s production 
processes g + g — > s + s and q + q — » s + s , in lowest order are well known: 




Figure 10. QGP strangeness production cross sections 0: Solid lines qq —> ss; 
dashed lines gg — » ss. a) for fixed a s = 0.6, m s = 200 MeV; b) for running 
a s (y/s) and m B (y^s), with a s (M z ) = 0.118. m s (M z ) = 90 MeV, m s (lGeV) ~ 
2.lm s (M z ) ~ 200MeV. 



Considering that the kinetic (momentum) equilibration is faster than the process 
of strangeness production (chemical equilibration), we use equilibrium Fermi-Dirac 
or Bose particle distributions f{pi,T) to obtain the thermally averaged strangeness 
production rate: 

, > _ J d 3 Pl J d 3 P2 a 12 v 12 f( Pl ,T)f(p 2 , T) 

{aVl ' cl)T = j<p Pl jd3 P2 m,TMp 2 ,T) ■ (5) 

This leads to the (Lorentz invariant) reaction rate per unit volume and time: 

A^ ss = \pl{t) {ov)^ sS , A"^' s = p q (t)p q (t){av)f- sS , (6) 

A ss^g 9 m = ps{t) ps{t) {av y T ^99m_ 

The factor 1/2, introduced here for two gluon processes, compensates the double- 
counting of identical particle pairs, arising since we are averaging considering all 
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reacting particles and thus each pair is counted twice. This invariant rate equation 
10 is the source term for the strangeness current: 

d^jg = ^ + = A "^ SS + A9 ^ SS ~ A sS -"^. (7) 

In the local rest frame of reference (v — 0) and ignoring the collective flow of 
matter we have: 



dt dt 2 Hg 



(t) (av)^ ss + p q (t)p,(t)(av)f^ - Ps (t) Ps (t) M« 9 *. (8) 



Evolution for s and s is identical, which allows to set p s (t) = ps{t). Using detailed 
balance to simplify, we obtain: 

dps = A (l- J$L} ) A = A"^ sS + A q ^ sS . (9) 



dt V P 2 s(°°) 

The generic solution at fixed T (p oc tanh(t/2r s )) implies that in all general cases 
there is an exponential approach to chemical equilibrium 

£fW _> i _ e -t/T. T =1 no) 

PT ' S 2yt(t)(av)^ ss + p q (t)p,(t)(av)^ ss + ...' 1 ' 

where the characteristic time constant r s is the ratio of the density we are 'chasing' 
with the rate at which the case occurs. There could be additional strangeness 
formation processes as is suggested in equation <|10[) in denominator for r s . 

The strangeness chemical relaxation time is rather short, evaluation of the diverse 
expressions presented here produces the result seen in figure 1111 where the large 
uncertainty is due to 20% uncertainty in the value of the strange quark mass. 
We recognize that the standard QCD cross sections <Jq^ s and particle densities 
associated with a temperature T > 200 MeV yield r s similar to lifespan of the plasma 
phase. Here, we note that both t s and the lifespan of plasma decrease with increasing 
temperature of the initial state: for the plasma lifespan this is due to more explosive 
outflow of matter driven by greater initial pressure, while for t s this behavior is result 
of both the increased particle density and increased average reaction energy. 

The temporal evolution of the strangeness density can be arrived at considering 
entropy conserving expansion T 3 V = Const, and allowing for the dilution of 
the strange quark density by collective matter flow. Introducing the phase space 
occupancy factor 7P GP and limiting ourselves to consider Boltzmann distribution, we 
find that 7^ GP evolves in time according to: 

s dt{ dT + T K 2 (z)J ^ 7s > ' ls {t) ~ „«. ' Z T 

where Ki are the Bessel functions. We see that the time dependence of temperature 
determines the scale at which 7p- GP evolves. Solutions of this equation show that 
one can oversaturate the final QGP phase space, since the large initial strangeness 
abundance, even if sub-equilibrium at that time, can exceed the equilibrium population 
at lower temperature |20| . 

The same methods can be applied in the study of thermal charm formation. The 
relaxation time constant seen in figure 1121 is large, indicating that charm will not 
be significantly produced in thermal parton collisions. On the other hand, charm 
is produced abundantly in first hard parton collisions. Benchmark values are 10 
cc pairs in central Au-Au at RHIC-200; and 200 cc pairs in central Pb-Pb LHC- 
6000 reactions 021- This yield is greater than the expected equilibrium yield at 
hadronization of QGP. 
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Figure 11. Strangeness thermal relaxation time constant as function of 
temperature. Dotted line is for fixed a s = 0.6 and m s = 200 MeV while the 
solid lines are for running QCD parameters, see figure [ToT bO. Thin solid line 
indicates the remaining uncertainty in a s while the hatched area indicates the 
domain for 20% variation of m s . 
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Figure 12. Charm thermal relaxation time constant obtained for a s (Mz) = 
0.118, m c (M z ) = 0.7 ± 7% GeV. 

4. Statistical hadronization 

Hadron production is well described in a large range of yields by the mechanism of 
statistical hadronization 03] ■ This approach works in high energy nuclear collisions 
far better than has been anticipated by Fermi 03], the inventor of this method. 
An important result of statistical hadronization is that within a particle 'family', 
particle yields with same valance quark content are in relation to each other thermally 
equilibrated, e.g., the relative yield of A(1230)(qqq) and N(qqq) or K*(sq) and K{sq) 
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Table 2. Four quarks: s, s, q, q — > four chemical parameters. 



7i controls overall abundance 
of quark (i = q, s) pairs 



Absolute chemical 
equilibrium 



A 4 controls difference between 

strange and non-strange quarks (i = q, s) 



Relative chemical 
equilibrium 



are completely controlled by the particle masses mi , statistical weights (degeneracy) 
gi and the hadronization temperature T. In the Boltzmann limit one has: 

n*_ = g*m* 2 K 2 (m*/T) 
11 g m 2 K2(m/T) 

Measurement of hadron resonances can thus sensitively test the statistical 
hadronization hypothesis. Since it is possible that the decay products of resonances 
rescatter, in principle not all resonances can be reconstructed by the invariant mass 
method. Thus, one expects that the experimental result must be below the ratio 
equation (|12fl expected for a given temperature. This will be the case for our study of 
K* /K and A(1520)/A for which the data is available. In fact, in our study of particle 
yields in section^ the experimental ratio gS] \M E3 : (K* + K*)/K~ = 0.26 ± 0.07. 
We find that this result is exactly in agreement with our hadronization temperature, 
T = 145 MeV, see table We also find A(1520)/A = 0.053 which is about twice as 
large as the experimental value |48|: 0.025 ± 0.07. This reduced experimental yield 
confirms that d-wave resonance A(1520) is particularly fragile in matter |49j. 

The yields of particles are aside of temperature also controlled by their fugacity 
Tj = e CTi / T , where <7j is particle T chemical potential. Since for each related particle 
and antiparticle pair we need two chemical potentials, it has become convenient to 
choose parameters such that we can control the difference and sum of these separately. 
For example for nucleons and antinucleons N, N the two chemical factors are chosen 
as: 

a N = n b + T In Jn, cr w = -fi b + Tin j N , (13) 

T N = lN e^/ T , T w = lN e-^' T . (14) 

The role of the two factors can be understood considering at the first law of 
thermodynamics : 

dE + P dV - T dS = a N dN + a w dN, 

= n b (dN - dN) +T In j N {dN + dN). 

The (baryo) chemical potential [ib, controls the baryon number, arising from the 
particle difference. 7, the phase space occupancy, regulates the number of nucleon- 
antinucleon pairs present. At the quark level (combining the light quarks u, d in one 
q, we have the situation shown in table |21 

There is considerable difference how the two types of chemical factors influence 
particle yield equilibration. This is best understood considering strangeness in the 
hadronic gas phase, the two principal chemical processes are seen in figure IT3l The 
redistribution of strangeness among (in this example) A, ir and N, K constitutes 
approach to the relative chemical equilibrium of these species. The production 
processes, on right in figure 1131 are responsible for absolute chemical equilibrium 
of strangeness. Achievement of the absolute equilibrium, 7 — ► 1, require more rarely 
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occurring truly inelastic collisions with creation of new particle pairs. These processes 
are absent in usual nonrelativistic chemical environments. 




Figure 13. Typical strangeness exchange (left) and production (right) reactions 
in the hadronic gas phase. 

When statistical hadronization of, e.g., a deconfined state of QGP occurs, the 
hadron yields are following closely the magnitude of the accessible phase space 
characterized by these chemical factors (and temperature). An important feature 
is the apparent under or over population of the resulting phase space. Namely, even if 
there is the pair abundance equilibrium in the primary QGP phase, in the secondary 
phase (here hadronic gas) phase space has different size in general. Absolute chemical 
equilibrium in the secondary phase requires a period of pair production accompanying 
hadronization. A long lasting QGP-HG mixed phase is today ruled out by experiment, 
and we should expect that the phase space occupancy chemical factors ii will be 
discontinuous at the phase transformation, and in general in the observed final state 

There are many different hadrons, and in principle, we should assign to each a 
chemical potential and than look for chemical reactions which relate these chemical 
potentials, e.g., on left in figure ITTA we infer /lA + zi^ = hn+Hk- However, more direct 
way to accomplish the same objective consists in characterizing each particle by the 
valance quark content |50|. forming a product of chemical factors, e.g., for p(uud), 

TCp(W) = itld X 2 u X d , Yp(uud) = 7«7<J A. u 2 A d , 

note here that: 

A; = e^ /T , fl g = -{n u + fi d ), X 2 q = X u X d X b = Xg. 

This implies relations between quark based = u,d, s and hadron based fj,i,i = b, S 
chemical potentials: 

q - 1 x - A « 

Mb — ^Mg Ms — gMb — MS, *s — 

An important anomaly is the historically negative S-strangeness in s-hadrons. e.g.: 

T A = luldls e^ + ^+^' T = 7u7d7,e ( ^-^ )/T , 
T x = luldls e { -^-^-^ )/T = i u i d i s e (-"»+w)/ r . 
The phase space density is: 

(15) 

-^-bosons <- gTOi/T ^g-j 
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and thus the 4-7T particle yield is proportional to the phase space integral, for example 
for 7r, N and JV: 

^tt ~ y d 3 P 1 2 . m,/T„ /1 fi) 2 

•^ = Cg V(2^ 7 -W"^-i ' lq ( } ' 

N _ f d 3 p 1 N _ f d 3 p 1 

77 - C 9n ttt-^— — 77 — L' 9n i 



V ""J (2tt)3 1 + j- 3 X- 3 e E / T ' V ™J (2^) 3 1 + 7g - 3 A+ 3 e s / T ' 

One can show that there is no influence of matter flow dynamics on this complete '47r' 
particle yield. 

When a small region of rapidity is considered, but the yield of particles is 
practically constant as function of rapidity, it is possible to imagine that the yield arise 
from a series of fireballs placed at different rapidities, and thus in this limit we also can 
proceed as if we had a full phase space integral. However, the proportionality constant 
C is, in this case in particular, but also more generally for a dynamically evolving 
system, not the system volume. We recall that though the individual phase space 
integrals are easily evaluated, to get the full yield one has to be sure to include all the 
hadronic particle decays feeding into the yield considered, e.g., the decay K* — > K + ir 
feeds into K and it yields. This actually constitutes a book keeping challenge in study 
of particle multiplicities, since decays are contributing at the 50% level to practically 
all particle yields, sometimes the decay contribution can be dominant, as is generally 
the case for the pion yield, yet each resonance contributes relatively little in the final 
count, it is the large number of resonances that matters. 

It is often more appropriate to study ratios of particle yields as these can be 
chosen such that certain physical features can be isolated. For example, just the two 
ratios 



Ra 



A + S + £ j = 021 = \~ 2 \- 4 = e 2Ms/T e -2 Mi) /T 



A + E° + S* H sqq 



-+- a* H ssq ~ q 

lead to a very good estimate of the baryochemical potential and strange chemical 
potential f3El , and thus to predictions of other particle ratios. The sensitivity to phase 
space occupancy factors ji derives from comparison of hadron yields with differing q, s 
quark content, e.g.: 

E-(dss) 7d7s 2 g~\ d \ 2 s S (dss) ldl 2 s gsX^Xj 2 



oc — — — — — (x 

A(dds) 7^ g A X d X s ' A{dds) l 2 d ls gAX^K 1 ' 

Note that i 2 q = i u i d and i u ~ 7d . 

An interesting application arises when we consider product of the above ratios, 
as than the result depends on is /id and temperature. In figure [21 we see how this 
product of ratios shows for the RHIC130 results that the equilibrium value ls/lq = 1 
(where we set i d ~ i u ~ i q ) is not compatible with the experimental results. Recall 
that the phase space occupancies we consider here are those arising from hadron yields 
and thus are applicable to the hadron yields, the excess above chemical equilibrium 
does not imply that the underlying state of, e.g., deconfined quark-gluon plasma is 
not chemically equilibrated. 

It is at first hard to believe that the overpopulation of the strange quark phase 
space, as shown in figure should be 50% grater than that of light quark phase 
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Figure 14. Evaluation of results of RHIC130 experimental particle ratio results 
showing that "fs/'Yq > 1 El- The cross is result of a global fit to all data, the 
bands reflect on the experimental uncertainty. Hatched area is the compatibility 
region for the three particle ratio products considered. 



space. However, this should have been expected! As the kinetic theory of strangeness 
production presented in section |3| shows, the lower mass quarks approach equilibrium 
yield faster. Thus, while practically all QGP produced strange quark pairs escape to 
be observed, light quark pair can be more easily reequilibrated in their abundance. 
This argument assumes that the hadronization is not entirely sudden. 

In the limit of very rapid hadronization, another mechanism favors abundance of 
strange quark pairs over light quark pairs in that there is a relatively small upper limit 
on 7q. As there is no time to expand the volume, hadron formation has to absorb 
the high entropy content of QGP which originates in broken color bonds. As is seen 
in figure Il5l the maximum entropy density S/V occurs for an oversaturated pion gas, 



lq 



a m„/2T 



1.6. Just below this Bose condensation condition, the entropy density 



is twice as large as that of chemically equilibrated gas. The entropy content of a 
non-equilibrium Bose gas has to be evaluated recalling 

' iPp d 3 x 



Sir 



{2irhY 



[(l + /,)ln(l + U)-U In U] 



1 



* e E w /T _ ! 



lq e 



i/m| +p 2 . 



(17) 



(18) 



We note that a similar constrain arising from kaon condensation limit yields the less 
restrictive condition: 



lq \\ 



±1 



< e 



(m K -m„)/T 



11 



(19) 
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5. Hadronization at SPS and RHIC 

Given yields of (strange) hadrons, we can apply the principles described above to study 
the hadronization conditions, and we would like to do this here as function of energy. 
Experimental results were available to consider RHIC Au-Au collisions obtained at 

y / s^" = 130 GeV and SPS Pb-Pb reactions at ,J~s%% = 8.75, 12.25, 17.2 GeV 

(projectile energy AO A GeV, 80 A GeV, 158 A GeV), as well as S-Pb/W 200 A GeV 

reactions occurring at \ fs$$ — 19.2 GeV. The procedure we adopt is to fit statistical 

parameters for Air particle multiplicity results. For SPS, the resulting % 2 /dof is 
shown in figure 1161 The solid squares are obtained allowing 5 parameters, i.e., 
full chemical noncquilibrium, the open triangles assuming that complete chemical 
equilibrium prevails. Open square results were obtained assuming that light quarks 
are in chemical equilibrium. Except at 40^4 GeV, we find in our approach involving 
chemical nonequilibrium yields and only statistical errors, statistically very significant 
fits. _ 

We use here SPS NA49 results [Hj, which include tt ± , K, K, A, A, 4> at 40, 80, 
160A GeV and at top SPS energy, we find that the experimental results for 3, 3, il, 
are reproduced nearly exactly even though for consistency we have not fitted these. 
Since this is a different data sample and we do not employ presently the central rapidity 
WA97 results, our fit result for some parameters will be slightly different, e.g., though 
we still report 7^ G > 1, the value we find is smaller than found using the complete 
data sample of NA49 and WA97. 
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Figure 16. Quality of statistical hadronization fit to hadron yields: \ 2 /dof for 
the 4 SPS reaction energies. 



Table 3. Results of fit of RHIC130 particle yields. Top line: statistical 
significance, following sections: fitted statistical parameters, hadro-chemical 
potentials, physical properties of the fireball. Stars indicate that a value is 
consequence of a constraint, and was not fitted. 



RHIC130 


non-equilibrium 


equilibrium 


X 2 /dof 


27/(21 - 3) 


230/(21 - 2) 


T 
\ q 

7 HG 

7! JG /7? G 


144.6 ± 1.3 
1.069 ± 0.008 
1.0198* 
1.62* 
1.53 ± 0.03 


169.1 ± 2.3 
1.067 ± 0.008 
1.0167* 
1* 
1* 


Mt [MeV] 
Ms [MeV] 


29.1 
6.9 


32.8 
8.1 


s/b 

E/b [GeV] 
S/b 

E/S [MeV] 


8.8 
32.8 
215 
153 


6.2 
31.8 
218 
146 


W s 


0.55 


0.49 



The RHIC fit favors non-equilibrium by such a large margin that we could not 
show the result on the same scale in figure El The actual values for x 2 /dof are seen 
in the top line of table |3] In obtaining these results it is assumed that 40% of weak 
decay cascades S — > A and A — > p are accepted whenever such corrections were not 
applied to data by the experimental groups. In the chemical nonequilibrium fit, we 
assume that the maximum value of 7 9 = e m *l 2T is attained, and thus this value is not 
fitted. The value of \ s is obtained from the requirement that strangeness balances, 
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(s — s) — 0. The other statistical parameter acquire values seen in table 01 We note 
that the hadronization temperature is greatly reduced in non-equilibrium approach 
compared to the value arrived at forcing chemical equilibrium. This is made possible 
by 7i > 1, in this case the high particle yields seen at RHIC can be arrived at a lower 
value of hadronization temperature. We also present the hadron gas parameters [is 
derived from X q , A s which can be used in back of envelope cross checks of expected 
particle ratios. These results constitute a slight update of our standard RHIC fit 5 , 
where the actual particle multiplicities are also given. 

Given the fitted set of statistical parameters (T,^j,7j) which characterize the 
phase space, we can evaluate, up to a common normalization constant, the strangeness 
(s) ~ (s), energy E, baryon number b, entropy S contained by all produced hadronic 
particles, as shown in table At the very bottom, we also present the Wroblewski 
ratio |51j . 

{da + uu) 

of newly produced strange quark pairs to all quark pairs. Strangeness production, 
while enhanced is not yet as abundantly produced as light flavors. We consider this 
ratio as function of y/s in figure 1171 The triangles for elementary p—p reactions are 
consistent with 'color string snapping' process which yields strangeness pairs with 
about 22% abundance of the light quark pairs, and nearly independent of y^i. In 
heavy ion collisions (circles), a much greater and \/s dependent result is seen. This 
indicates that a new strangeness production process arises, which could well be the 
thermal production of strangeness in QGP. We also see, in the background, the results 
of a similar analysis obtained within the chemical semi-equilibrium approach |52j . 
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Figure 17. Wroblewski ratio 1511 : only newly made s- and q-pairs are counted 
in comparing strange to light quark pair production. 



The nonequilibrium parameters "f s /jq, 7q are displayed in figure ITH1 The results 
for J s /jq, on top, show that only at RHIC j s > 7 9 . Open squares show that when 
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Figure 18. Fits of, top, 7s/7q and, bottom, 7 q chemical nonequilibrium 
parameters as function of collision energy. Open squares results for j q = 1, 
in which case the left hand side gives j s . 



we fix j q = 1 the SPS fits, in general, yield 7 S < 1, as is widely reported. However, 
bottom, we see that 7 9 > 1 in all cases, reflecting on the excess in charged hadron 
multiplicity. We note that the 40^4 GeV SPS result (lowest energy point) deviates 
from the behavior systematics of the other results. 

Since we know the incoming energy and we can evaluate the amount of thermal 
energy found in the final state hadrons, we have a good measure of the energy stopping 
for the different collision systems — to be precise in estimate of the energy stopping, 
we would have also to account for the energy transfered from intrinsic degrees of 
freedom to the collective flow. The behavior of the ratio of thermal energy to collision 
energy is seen in figure 1191 The result is in so far surprising as energy stopping is 
somewhat greater at RHIC130 than at top SPS energy. On the other hand, the rise 
of the energy stopping power towards lower SPS energies is expected, as is the steep 
rise for the asymmetric collision system S-Pb/W. 

With the relatively great variability of stopping power seen in figure H^l it seems 
more prudent to use in study of energy dependence of different physical properties 
of the hot fireball not the collision energy but the final state thermal fireball energy 
content. This is strongly supported by the study of the systematics of strangeness 
production per baryon, seen in figure |2"U1 On top, we show the strangeness excitation 
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Figure 19. Fraction of energy stopping at SPS and RHIC: results are shown for 
40, 80, 158A GeV Pb-Pb, 200A GeV S-W/Pb reactions and at RHIC for 65+65A 
GeV Au-Au interactions. 




function as usually shown, using ^snn as variable, and normalizing to the baryon 
yield. The baryon number is conserved in hadronization, and this result is directly 
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Figure 21. Entropy production per thermal participating thermal baryon S/b 
shown as function of the intrinsic thermal energy per baryon, E* \ N . 



telling us how much strangeness was available prior to hadronization. We recall 
that the true 'thermalized' participants are about 5-10% fewer than is inferred from 
geometric consideration. On bottom in figure I2TJ1 we present the same result using as 
variable the intrinsic energy content. There is a significant simplification of the result 
and we find 



Importantly, there is no sign of a deviation of the systematics of strangeness 
production. 

Similarly, considering the specific entropy content in figure |^ we find that 
entropy production per unit of available energy is universal: 



The entropy production mechanism which has this simple result for a wide range of 
intrinsic fireball energy is not known, but appears to be common for both SPS and 
RHIC energy range. 

Comparing the results seen in figure l2lH and figure l2"Tl we note that, as function of 
the intrinsically available thermal energy, the strangeness production rises somewhat 
faster than entropy, which expresses the fact that strangeness is more enhanced than 
hadron multiplicity comparing RHIC to SPS. Within QGP phase, the strangeness 
to entropy ratio s/ S, shown in figure 1221 characterizes the approach to chemical 
equilibrium of strangeness, as compared to saturation of the availability of quarks 
and gluons. Both quantities may increase slightly in hadronization, but are expected 
to be practically conserved. Since the rise of s/S from SPS to RHIC is relatively slow, 
we show the behavior as function of both, top ^/snn, and bottom of Ef^ N . There is 
clear increase of this ratio comparing SPS to RHIC, but the magnitude of the effect 
is somewhat dependent on the method of evaluation of chemical conditions. 

There is perhaps an indication, in figure 1221 of a possible increase of strangeness 
yield per entropy at the lowest energy, consistent with the somewhat anomalous 
rise of K + /ir + ratio reported by the NA49 experiment. Dynamically, this may be 
consequence of a longer lived QGP phase, which lacks the strength to undergo a 
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Figure 22. Strangeness per entropy: top as function of ^/satjv, bottom as 



function of . 



fast collective flow. Longer lifespan even at lower initial temperature would help 
strangeness production. More systematic study of strangeness production in this 
energy range could thus offer significant insights into the properties of the deconfined 
phase. Ultimately, as the energy is reduced below the QGP threshold, there should 
be a significant cut in strangeness yield. 

6. Onset of quark-gluon plasma formation 

The properties of hot QCD matter can be well described in terms of a QGP 
model which agrees with latest lattice results [^J 123 123- Thus, in principle, 
we understand the behavior of the QGP and can analyze the meaning of the relatively 
low hadronization temperatures. In figure l2"51 we show with solid squares the SPS and 
open square the RHIC hadronization points we obtained, dashed line corresponds to 
the cross over from the deconfined to confined phase. In the background, we see as 
triangles hadronization analysis for AGS and SIS assuming equilibrium conditions. 

Our results are tabulated in tabled obtained enforcing strangeness conservation 
(except for S-W) and allowing chemical nonequilibrium. As we go to press a fit at 
30 GeV is still not possible, with particle yields not available. Statistical errors on 
the fit results are small and are not shown, in figure 1231 we indicate our estimate 
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Figure 23. Hadronization conditions in T, fn, plane: RHIC130 open square, SPS 
solid squares. 



Table 4. The results of nonequilibrium fits at SPS and RHIC: for each collision 
system we present the collision energy in laboratory, reaction energy ^/sjvjVi the 
chemical freeze-out temperature T, and baryochemical potential fj,i,; these results 
are shown in figure l23l The bottom line gives the experimental ratio K/71-, where 
available, see text for details. 





Au-Au 


S-W 


Pb-Pb 


Pb-Pb 


Pb-Pb 


Pb-Pb 


V^'lab L«eVj 


65+65 


200 


158 


80 


40 


30 


v /s/Viv [GeV] 


130 


19.4 


17.2 


12.3 


8.75 


7.60 


T [MeV] 


144.6 


143.5 


148.5 


150.3 


139.5 




Mft [MeV] 


30.1 


165 


236 


302 


370 




K/tt 


0.16 




0.12 


0.11 


0.11 


0.12 



of the systematic SPS result error based on 5% individual data point error. The 
relative changes between the fit results to NA49 40,80,158 GeV reactions are most 
likely significant, since the systematic error is common to these results. Thus it is the 
80 GeV SPS hadronization point which is closest to the equilibrium phase transition 
between HG and QGP. The RHIC point (open square) though at low baryochemical 
potential could be at lower temperature than are the 160 GeV SPS and 80 GeV SPS 
results. The lowest temperature is found for the S-W point, but this result is obtained 
within a very different procedure, using central rapidity results of WA85 and a flow 
model. If at all, the temperature of hadronization seems to rise as the collision energy 
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falls, until we nearly touch the equilibrium phase transition condition. From thereon, 
we record a different behavior, with hadronization condition at 40 GeV moving more 
inward the hadron domain, with AGS and SPS equilibrium points deep within hadron 
gas phase. 

To better understand this result, we must remember that the fireball is not a 
piece of deconfined matter sitting still. It undergoes a rather rapid collective explosive 
flow. Collective motion of color charged quarks and gluons contributes an important 
collective component in the pressure, as can be seen considering the stress portion of 
the energy-momentum tensor: 

T«=P% + (P + E )J^L-. (21) 

The rate of momentum flow vector V at the surface of the fireball is obtained from 
the energy-stress tensor Tkf. 

p = f-n = Pn+(P + e)^^. (22) 

The pressure and energy comprise particle and the vacuum properties: P = P p — B , 
e = e p + B . 

The condition V = reads: 

Bn = P p n+{P p +e p ) V -f^. (23) 
Multiplying with ft , we find: 

B = P p + (P p +e p )-^L, «=^#. (24) 

i - vi vi 

This requires P p < B: QGP phase pressure P must be negative when the dynamical 
pressure comprising effect of the flow runs its course, V — > 0. A fireball surface 
region which reaches V — ► but for which v ^ 0, i.e., it continues to flow outward, 
is torn apart in a rapid filamentation instability. This situation can only arise since 
the quark-gluon matter presses again the collective vacuum which is not subject to 
collective dynamics. 

Phase boundary between the hadron gas domain and quark-gluon plasma, which 
includes the effect of the 'wind' of flow of QCD matter is shown in figure El for a 
geometric model with k — 0.6. With increasing flow velocity the phase boundary 
set at V — > moves to lower temperatures. This effect is larger near to fib = 
than the effect fit, has on change in T. We can now understand the behavior seen in 
figure at RHIC there is the largest supercooling and the flow velocity is greatest, 
yielding smallest hadronization temperature. As collective flow velocity decreases, the 
hadronization temperature increases, though this effect is somewhat compensated by 
the influence of the the increase in fit which implies a decrease in T. The impact 
of the flow velocity is smallest near to 80A collision energy, where the hadronization 
temperature is nearly reaching the HG-QGP phase boundary. 

This would not suggesting that the onset of QGP formation is in the collision 
energy window 40-80A GeV. This would be in conflict with the result shown in last 
line of table where we present the experimental result for 



K/tt = -y/K + /7r+ x K~/tt-. (25) 

The advantage of taking the geometric mean of the charged kaon to pion ratios is 
that the chemical potentials of strangeness and baryon number cancel in Boltzmann 
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Figure 24. Hadronization boundary in the /Jj,,T plane: solid (point hadrons) 
and dashed (finite volume hadrons) lines are estimates of QGP-HG boundary for 
a system at rest. Dotted lines are for finite size hadrons and include collective 
flow velocity for v 2 = 0, 1/10, 1/6, 1/5, 1/4, 1/3 . Thick solid line: breakup with 
v = 0.54, k = 0.6 1411 . The measured point relates to the hadronization condition 
at 158A GeV. 

approximation. The constant value we see at SPS leaves no space for any discontinuity 
in strangeness production in SPS energy range, including the newly available result 
at 30 GeV ^21- This finding contradicts the claim that the rise of K + /tt + ratio with 
decreasing collision energy suggests new physics occurring near 30A GeV [S7| (but 
does not contradict the fact, which is claimed on false grounds). The consideration 
of K/7r indicates that the variation of the baryochemical and strangeness chemical 
potentials fully account for peak in K + /ir + production. 

The increase in K/tt ratio between RHIC and SPS energies is an indicator of 
new physics, and an energy scan between RHIC and SPS energies would at first seem 
appropriate. However, the 40% increase we see moving on from SPS to RHIC has 
been in this paper understood as the result of the initial state changes and faster, 
more explosive expansion of the fireball. Still, it would be good to make sure that this 
rise is not step-like, which would suggest a phase change. 

Our findings make the issue where is the onset of QGP formation at first 
enigmatic. We have seen that specific per baryon strangeness and entropy are both 
evolving smoothly in the energy range beginning at AOA GeV and through the RHIC 
energies. We see in figure |231 that near 80 A GeV supercooling seizes, and below this 
collision energy the particle freeze-out occurs in HG phase This means that if QGP 
phase is formed below 80A GeV, there is reequilibration of particle yields. Still, we 
expect that strangeness yield would be preserved, and any enhancement of strangeness 
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to entropy yield as expressed by the K/V ratio equation (|25|l should remain visible. 
The top AGS energy at 11. 6 A GeV for the Au beam yields [SHI: 

K/tt = 0.077. 

An increase of AGS result by 55% is required to reach the SPS level present at 
30A GeV. This increase is sufficiently sudden to suggest that a change in reaction 
mechanism must occur between 30 and 10A GeV beam energy. Since the only change 
we are expecting is the onset of deconfinement, we have all reason to hope that a 
rather sudden rise in K/-7T can be observed in this energy domain. Naturally we 
also expect that evolution of specific strangeness and entropy production will show a 
sudden increase in this energy domain. 

7. Anything new at LHC? 

An often posed question is if there can be a future for strangeness as signature of 
QGP at the LHC energy scale which is 30 times greater than at RHIC. In principle, 
one would think that the extreme conditions expected will reduce the importance 
of strangeness as diagnostic tool, since the difference between u, d and s flavors will 
diminish. We do not believe that charm and bottom flavors will 'replace' strangeness 
as an observable, rather these new degrees of freedom produced in first collisions will 
complement the physics potential of strangeness, which as we now argue, will further 
grow at LHC. 

The QGP based production processes at RHIC lead to a highly oversaturated 
strangeness phase space yield. The mechanisms responsible for this are likely to 
augment as the energy available in the nuclear collision increases: it is probable that 
at LHC we reach greater initial temperatures and more explosive transverse flow. As 
the system volume expands, but entropy remains nearly conserved, there is significant 
reduction of particle pairs. Since the reannihilation of light quarks is favored by 
cross section, compared to the reannihilation of more massive strange quarks, this 
faster expansion is driven predominantly by the energy derived from light flavors. 
Strange quarks, once produced should remain more if not mostly preserved till the 
hadronization. 

As a consequence, we expect that "f s /lq will further increase at LHC compared 
to the large RHIC value, which at yfs — 130^4 GeV is about 1.5 . An increase in this 
ratio triggers an increase in the yield of kaons, as compared to pions. At the nearly 
baryon free conditions prevailing at LHC, we can consider X q = X s — 1, hence the 
value of K ± /7r ± ratio in a statistical hadronization model is proportional solely to 
7 S /jq and depends on the hadronization temperature T. When 7,5 = 7 9 = 1, this 
ratio as function of T is shown in figure l2"B1 as dashed line. This chemical equilibrium 
result applies also at finite baryon density (SPS) to the geometric mean of the charged 
particle ratios, equation l|25|l . At T ~ 160 MeV, we see in figure [35] K/tt ~ 0.133. 

The solid line, in figure 1551 shows the maximum possible value of K/tt ratio, with 
both 7 S and "f q increasing to the maximum allowed by the requirement that neither 
kaon nor pions condensation arises, compare equation (|19|) : 

7s7 9 = e" K/T , J q l g = e m ~ /T , for X q = X s = 1. 

The arrow indicates where the lines will meet at very large hadronization temperatures. 
The exact location of this limit is in part result of yet limited knowledge of hadronic 
high mass resonances. In fact, it can be suspected that the slight decrease in the 
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Figure 25. Kaon to pion ratio (see text) as function of hadronization 
temperature for chemical equilibrium (dashed line) and with maximum allowable 
7„ and 7 9 (solid line). 

ratio K/7r with increasing hadronization temperature is spurious, originating from our 
lack of knowledge of the hadron strange and non-strange mass spectra. Moreover, it 
can be expected that hadronization temperature remains at or below the Hagedorn 
temperature T ~ 160 MeV at LHC. 

The important result we see is that K/ir can almost triple without violating any 
fundamental principles, in a scenario in which initial thermal and very hot QGP phase 
is formed and expands explosively. Of course, this does not prove that this will happen 
at LHC, though a significant increase of K/tt should occur. Perhaps more spectacular 
is the expected rise in the Wroblewski ratio W s , which compares the effectiveness of 
strange quark production to light quark production |51j . counting only newly made 
quark pairs, compare equation (|2*U|) and figure IT7I W s ~ 0.2-0.25 in p-p reactions 
expressing the low yield of strange quark pairs available. In heavy ion collisions, values 
as large as triple this result have been observed, expressing, as discussed, the large 
enhancement of strangeness yield as compared to light quark yield. Somewhat smaller 
W s values are expected in chemical equilibrium in baryon-free matter, as shown by 
dashed line in figure The solid line, in figure l2lfl shows the great enhancement of 
W s should maximal overpopulation of strange quark phase space indeed be established 
(note logarithmic scale in figure l2*6l . The arrow indicates the large T limit of W s which 
against naive expectation is found somewhat below unity. This is an expression of the 
asymmetry in the number of undiscovered strange and non-strange high mass hadron 
resonances. Again, this is a purely academic point, as hadronization is expected at or 
below the Hagedorn temperature T ~ 160 MeV. 

The new physics presented in this report, include a consistent chemical non- 
equilibrium analysis of the SPS and RHIC experimental results which shows great 
consistency of physical properties and reaction mechanism in these two domains. We 
have demonstrated that the great enhancement of strangeness production at RHIC 
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Figure 26. Wroblewski ratio W s for chemical equilibrium (dashed line) and 
maximum allowable 7„ and ■yq in baryon-free matter. 



is originating in the same mechanism as strangeness production at SPS. We have 
illustrated how low hadronization temperatures derive from supercooling due to the 
fast collective flow of deconfined matter, and we have above presented a possible 
scenario for a large strangeness anomaly at LHC. 
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